REDUCED FORMULATION OF A STEADY FLUID-STRUCTURE 
INTERACTION PROBLEM WITH PARAMETRIC COUPLING 

TONI LASSILA* AND GIANLUIGI ROZZAt 

Abstract. We propose a two-fold approach to model reduction of fluid-structure interaction. The state 
equations for the fluid are solved with reduced basis methods. These are model reduction methods for 
parametric partial differential equations using well-chosen snapshot solutions in order to build a set of 
global basis functions. The other reduction is in terms of the geometric complexity of the moving fluid- 
structure interface. We use free-form deformations to parameterize the perturbation of the flow channel 
^^1 at rest configuration. As a computational example we consider a steady fluid-structure interaction 

Cw ' problem: an incmpressible Stokes flow in a channel that has a flexible wall. 
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1. INTRODUCTION 



Many problems in fields such as aerodynamics and biomeclianies can be expressed as fluid-structure 
interaction (FSI) problems. The mathematical modelling of such coupled problems consists of four 
main parts: solution of the fluid equations given the current fluid geometry, solution of the structural 
r^ [ displacements given the normal stresses exerted by the fluid, fulfillment of coupling constraints to achieve 

force balance across the interface, and transport of the fluid-structure interface. Even when both the fluid 
and structure equations are independently linear, the geometric variability of the fluid-structure interface 
yields a significant nonlinearity in the coupled system. Solution methods for fluid-structure interaction 
problems are iterative in nature and involve the repeated solution of the fluid and structure equations in 
many different configurations. 
^ I Fluid-structure interaction problems arise in modelling of the arterial deformations in the human 

cardiovascular system. A comprehensive treatment of FSI strategies in cardiovascular modelling can be 

P^ I found in [13l [15] . Blood flow is pulsatile and the displacements of the arterial walls are relatively large. 

fT^ ■ This leads to stability considerations that necessitate the use of implicit solver strategies of the coupled 

FS problem, again increasing computational complexity. 

(^ • Model Order Reduction (MOR) is a cross-disciplinary field that strives to systematically reduce the 

f^ I dynamics of a complicated ordinary or partial differential equation model to a simpler and computation- 

ally more tractable one. Model reduction techniques have also been proposed for nonlinear systems. One 
particularly common approach is the Proper Orthogonal Decomposition (POD) technique that aims at 
decomposing the dynamics of a time-dependent system into fundamental modes, and then choosing only 
. ^ very few of the most important modes to represent the entire dynamics of the system, 

j^ ■ A modal reduction of the structural equations was adopted in |llj to reduce the complexity of an FSI 

problem involving pulsatile flow in a channel, but no attempt to reduce the complexity of solving the fluid 
equations was made. In another work by the same author |10j the steady problem was parameterized using 
shape functions for the boundary trace of the pressure obtained from an assumption of near-Poiseuille 
flow proflle. The strong coupling of normal stresses on the interface was treated in the reduced space 
obtained by considering only the leading structural eigenmodes, and solved as a least-squares optimization 
problem. 

Another approach consists of the Reduced Basis (RB) methods, which is generally used on parametrized 
partial (and ordinary) differential equations. RB methods are based on a greedy sampling algorithm that 
chooses snapshot solutions of the parametric PDF at different parameter values for constructing a global 
approximation basis. For an introduction, see [M] and [20]. The reduction of steady incompressible 
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Figure 1 . Stokes flow in a channel with a flexible wall section 

Stokes flow (without FSI) in an arterial bypass configuration by parameterization with piecewise affine 
maps was considered in [21] . The same problem parameterized with nonaffine maps was considered in [19) 
for Stokes equations and in [16] for Navier-Stokes equations. We also mention the reduced basis element 
method proposed for Stokes flows in [8] that was derived to address complex flow networks with multiple 
branches and junctions. A posteriori error estimation for Navier-Stokes equations solved by reduced basis 
methods is considered in [31 [53] . We use reduced basis methods to approximate the solution of the Stokes 
equations in a parametric flow geometry. 



2. PROBLEM OF STEADY STOKES FLOW IN A FLEXIBLE CHANNEL 

We denote by H™{X) the usual Sobolev space of real- valued functions on X with mth weak derivatives 
being square-integrable. Let Oq C M^ be a bounded domain with Lipschitz boundary dflQ that represents 
the fluid domain at rest configuration (without deformations induced by the fluid flow). We denote the 
flexible part of the boundary as Eq C dflo- We assume the geometry displayed in Fig. [T]- a straight 
2-d tube, where the displacement of the upper wall is described with a function i] G 7Jg(Eo) := {4> G 
H^CEq) : (j){a) = 4>{b) — 0}. The tube is at "rest" configuration when 77 = 0. We denote by 0(77) the 
deformed domain. In the deformed domain we have the incompressible steady Stokes problem to find 
u e [H\n{f]))]^ and p e Ll{n{r])) s.t. 



(1) 
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on I](r?) 



where ly is the viscosity, and / is the volume force acting on the fluid. The Dirichlet data are assumed 
to satisfy a conservation principle 



(2) 



uq ■ n dT = 0. 



dn{ri) 



The fluid equations are then to be coupled together with the structural equations given the coupling 
condition that the displacement of the structure is the same as would be obtained from the structural 
equations, where the RHS is the stress exerted by the fluid on the interface. The structural equations 
can be chosen in several different ways (in "W and fVD different fourth order models are considered). We 
choose, after a suitable assumption of small structural displacement, to model the displacement of the 
wall using a 1-d elliptic equation ("elastic membrane") to find 77 e HqCSq) s.t. 



(3) 



/ K{x)ri'(t,' dV^ j T{u,p)(l)dT, V^eiJcJ(So) 



where r(tt,p) is the normal traction on the interface imposed by the fluid 



(4) 



t{u,p) = [pn — V (Vm + Vu ) • n\ 
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obtained from the solution {u,p) of the fluid equation ^, and K{x) > Kq > is the spring constant of 
the waU to displacements in the normal direction. 

In operator form, the coupled fluid-structure problem can be considered to be of the form 



(5) 



Because the fluid domain depends on the domain defined by the function 77, this is in fact a nonlinear 
free-boundary problem. 
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3. PARAMETERIZATION WITH FREE-FORM DEFORMATIONS 

To reduce the complexity of the free-boundary problem ([5|) we introduce a parameterization of the 
fluid domain, ri(/x), where the parameters /x G I? belong to some low-dimensional parameter space. One 
option to construct Q{n) would be to use a family of parametric curves, such as B-splines, to directly 
parameterize the displacement function 77. We prefer to use free-form deformations [22] . 

Free-form deformations are a technique of obtaining small, parametric deformations of an arbitrary 
reference domain flo C M"^ for d = 2,3. First we map the reference domain to a subset of the unit 
square with a continuous invertible map $ : f2o — > (0, 1) x (0, 1) and then overlay on ^(JIq) a regular 

The 



grid of control points P^ 



0, . . . , L and m — 0, . . . , M, so that P^ 



[i/L m/M] - 



parameterization is obtained by allowing a subset of the control points to move and then using a spline 
basis to construct a smooth deformation map as a function of the positions of the control points. The 
perturbed control points are given by a set of {L + 1){M + 1) parameter vectors /i^^„j. For each fi, a 
parametric domain map T : JIq x 2? — > ^(m) is defined as 



(6) 



\e=0 m=Q J 



iL,M 



where hg ^ are suitable polynomial basis functions. Here we take 
that is the tensor products of the 1-d Bernstein basis polynomials 

bfis) = (fi (1 - .)(^-^).^ b'^it) = (^ (1 - tr^-h- 

defined on the unit square with local variables (s,i) G (0, 1) x (0, 1). It is typical that only some of the 
control points are allowed to move freely, we do not use the full parameterization but only a few of the 
most relevant shape parameters. To parameterize the flow channel we use a 10 x 2 grid of control points. 
Six control points in the top row are moving in the a;2-direction, resulting in a smooth parameterized 
deformation of the channel with 6 real parameters. The result when applied on a sample computational 
mesh is shown in Fig. [2] 

After the parametric domain fl{fi) is obtained as the image of a fixed reference domain flo under the 
map T : i7o X I? — ^ ^(m); which is assumed to be differentiable and invertible, we proceed to transform 
the PDE on a parametric domain to a parametric PDE on the fixed reference domain. Denoting by 



VxTri the Jacobian matrix of T^ we define the transformation tensors for the viscous term 



(7) 



VT„ 



Jt^ Jt^ A^^{Jt,) 



and the pressure-divergence term 

(8) 



XTr, -.^ Jt^^ detiJrJ 



respectively. The inhomogeneous Dirichlet conditions are handled by lifting to a space of proper boundary 
conditions the solution of the homogeneous Stokes problem on the reference domain to find u(ri) G 
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Figure 2. Free- form deformation of the flexible channel with six parameters 

dUkr , 

dxi ^ dx^ 

duk 



(9) 



duk. dvk r , dvk\ 

"T^^^T.kjj^+PiXTjkj-g^j dno - ^^ 



detiJrJif^ + fmt]k dno, 



/ 1 [XTjfcj 



dxi 



dno = 0, 



yqeLliQo) 
with summation understood over the indices i — 1,2 and j — 1,2. 

4. PARAMETRIC FLUID-STRUCTURE COUPLING ALGORITHM 

We can alternatively formulate the fluid-structure problem ([5]) as a calculus of variations problem to 
find 1] G Hq{^q) the minimizer of 



(10) 



mm 



J So 



s.t. 



-lyAu + ^p = f in VL{ij) 

V • M == in 17(77) 

If a solution of ([5]) exists then it is also a minimizer of (|10p . Problem ([TU)) in the parametric form given 
by © is still valid: find fi E V s.t. 



(11) 



mm 



\Kd,^7j + Tiu,p)\^dT, 



s.t. 



-vAu + Vp= f in ri(/x) 

V • M = in ri(/x) 

but this time we expect that the optimal value of the functional is J* > and the coupling of stresses 
on the interface between fluid and structure is only achieved in a least-squares sense. The "goodness of 
fit" depends on the dimension of the parameter space T). We call this the strong parametric coupling 
approach. 

It should be noted that we perform a reduction by parameterizing the displacement 77 of the fiuid- 
structure interface. In ^10) the traction r (which corresponds to the boundary trace of the pressure in 
this case) was instead parameterized. This alternative approach allowed the author to solve directly a 
parametric minimization problem of the form (jlip . Instead, we must consider the minimization problem 



(12) 



mm 



|r;(/x) - 77|2 + |9,rK/x) - 5,f)pl dP, 
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S.t. 



-vAu + ^p= f^ in rj(/i) 



V ■ M = in fi(/x) . 

^Kdxr,f]^T{u,p)^() on E(/x) 

The reason for this is that, in order to reduce the complexity of the fluid equations, we must be able to 
transfer the parametric dependence to the Stokes equations. This can be done more easily by parame- 
terizing the displacement rather than the boundary trace of the pressure. 

A computational algorithm for the solution of ([T2l) is detailed in Algorithm 1. This is a fixed point 

Algorithm 1 Strong parametric coupling of fluid and structure 

Require: initial guess /x*^ 
Let A: = 0. 
repeat 

Solve the discretized fluid equations for {u'^,p'^) in n{fjb'^). 

Compute assumed traction f = p'^ls(^) — i' [Vit'"' + [Vtt'^]'^] from [u'^.p'^). 

Solve the minimization problem 

/x'^+i := argmin \ f {|7,(/x) - fif + |a,77(/x) - d^fif} dT 

to obtain a new configuration parameter, where ri{fi) is the displacement of the interface given by 
the parameterization of the geometry and the assumed displacement i) £ ffp (Eq) is the solution of 

/ Kf]'(j)' dV ^ I T(j)dT V e Ho (So) 

6: Set k := k + 1. 

7: until stopping criteria |/x'^+^ ^ A*'^| < e is met 

algorithm in the parameter space. The algorithm's computationally most intensive part is the solution of 
the parameterized fluid equations. We next discuss the reduced basis method for approximation of the 
fluid solution {u^ ,p^). 

5. REDUCED BASIS APPROXIMATION OF THE STOKES PROBLEM 

The computational efficiency of solving the coupling problem ()12|) hinges on the efficient solution of 
the fluid equations on the parametric domain r2(/x). We present briefly the general approach to reduced 
basis approximations for the steady Stokes problem. For more details see [12] and pi] . 

We denote by y := H^{VLo) x H^{^o) the velocity space, Q := L^(rio) the pressure space, and the 
continuous parametric bilinear forms ^:3^x3^x2?->R and B : Q x y x 1) ^ M.. The weak form 
(|9]) of the parametric Stokes equations for velocity and pressure can be written: given any /x £ 2?, flnd 
m(/x) £ y and p(/x) £ Q s.t. 

A[u, v; /x) + B{p, v; /x) = {F{p),v)y for all t; £ 3^ 

^^^^ B{q, m(/x)) - (G(/x), q)Q for all q £ Q, 

where F{fi) ^ y' x V and G(/x) £ Q' x 2? are continuous parametric linear forms. The inhomogeneous 
Dirichlet condition is treated by introducing an extension of the boundary data tto £ H^{ilo) and look- 
ing for a solution u £ HQ{il.o) x HQ{no) x T> to p^ with a modified right-hand side {F{fi),v)y :— 
{F{fi), v)y — {A{fi)uo, v)y. The solution is then recovered as m(/x) = u{fj,) + Uq. 
Define a(/x) the parameter-dependent coercivity constant of the problem as 

(14) a(/x):=inf -^^"'"''^^ 



vey 



2 

\y 



and /3(/x) the parameter-dependent inf-sup constant of the problem as 
(15) /3(/x):=inf sup ^lf';;;^^ 
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A necessary condition known as the Ladyzhenskaya-Babuska-Brezzi condition for problem (J13l) to have a 
unique solution is that a(/x) > ao > and /3(/x) > /3o > 0. We assume in what follows that the coercivity 
constant is always positive and concentrate on the inf-sup constant. 

From here on we forego any interest in the infinite-dimensional problem and assume instead that both 
y and Q are replaced by strictly finite-dimensional spaces yh and Qh- More concretely, we consider only 
the finite-dimensional problem obtained by discretization of the partial differential equations using the 
finite element method. The discretized model problem 

A{uh,Vh;fJ.)+B{ph,Vh;n) = {Fh{tJ.),Vh)y^ for aWvh eyn 

B{qh,Uh;n) = {Gh{fJ'),qh)QH for aU qh G Qh, 

involves a similar inf-sup condition to guarantee the existence of a unique velocity-pressure pair satisfying 



/-,-x a f \ ■ f B{qh,Vh;tJ') ^ « ^ n 

(17) P/i(m) = mi sup Y — 7—r — 7— >Po>0. 

gheQhVheyh \m\\Q\\Vh\\y 

In the case of the steady Stokes problem we employ the P^ — P^ Taylor-Hood elements [4j [17] that satisfy 
the condition (|17l) without the necessity of adding extra stabilizing terms. 

The reduced basis method seeks an approximation of the finite element solution to (|16p . The first 
step is to construct a set of global basis functions. Let fi^, . . . , fi^ be a set of snapshot parameter values 
chosen according to some rule. We denote by «(/.»") and p(/.t") the corresponding snapshot solutions 
for the velocity and pressure. The snapshot solutions are obtained as solutions of a stable finite element 
formulation of (TTSl) . The reduced basis velocity space is then defined as 



(18) span(u(/xi), . . . , m(/x~)) =: yN C yh, 
and the reduced basis pressure space as 

(19) span(p(/xi), . . . ,p{fi^)) =: Qn C Qh- 

We define the reduced basis Galerkin problem: for any /x G 2? find un G J^tv and pN G Qn such that 
A{uN{n),v;ti)+B{pN{fJ'),v;fi)^{Fh{fJ'),v) for ah u G J^at 

B{q,UN{n);n) = {GhifJ-),q) for ah q G Qn- 

Then (Mjv(/i),p7v(At)) is a reduced basis approximation for {uh{n),PhitJ'j) the solution of (fT6|) . To 
represent and solve these equations in matrix form, an orthonormal basis is constructed for yN and Qn 
to guarantee algebraic stability as iV — >■ cxd [?!] . 

Provided that the dimension N is chosen to be much smaller than the dimension of the finite element 
spaces, this system is inexpensive to solve. If the dependence of the bilinear forms on the parameters is 
smooth and the snapshot parameter values fi^ are properly chosen, we find that choosing a moderately 
small N permits rapid convergence towards the finite element solution. For details on how to choose the 
snapshots with the help of a greedy algorithm that controls the residual a posteriori approximation error, 
we refer to existing works on reduced basis methods [6l [141 [2Q] ■ 

A key issue in the reduced basis approximation of the Stokes problem is approximation stability. We 
remark that use of inf-sup-stable finite elements for the snapshot computation does not necessarily result 
in a stable reduced basis approximation when the basis functions are orthonormalized, as was observed 
in [3T]. On the other hand, the orthonormalization is vital for the algebraic stability of the reduced basis 
approximation. Thus there is a need to extend the concept of inf-sup constants into the reduced basis 
formulation. This is typically done by defining the so-called supremizers in velocity space and enriching 
the reduced basis velocity space yN accordingly. Following |21) we define the inner supremizer operator 
T*^ : Qh —>■ yh as the solution of the elliptic problem 

(21) {T''q,v)y^Biq,v;fj,) for all t;,, G iV,,. 

After the pressure approximation space Qn and velocity approximation space yN have been constructed 
we add all the velocity supremizers, i.e. velocity fields of the form T'^{p{fi"')), to the velocity space: 

(22) yN -.^yN U span[r'^(p(M^)),...,T'^(p(M^))] , 
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where y^ is a stabilized reduced velocity space. Then 

(23) /3w(m) :- inf sup ..^^'^''''^^ > Mfi) > /3o > 0. 

As a consequence of this approach to stabilization the dimension of the velocity approximation space grows 
to 2N, where N is the original dimension of velocity and pressure approximation spaces respectively. 

The computational benefits of the reduced basis method rely on the assumption that the parametric 
bilinear forms are afEnely parameterized: 

Ma Mt 

(24) A{u, v;,,)=J2 er (m)^'"(«, v), B{p, v; ,,) = Y. 0r(M)6'"(P, «) 

m—l m—1 

for some parameter- independent continuous bilinear forms A"^ and S'", and parametric coefhcient func- 
tions 8™ and 0™ that can be readily evaluated. Then the finite element equations (|16p split into 

M„ Mb 

Y^ er(/x)^"(wjv(/x),t;) + Y. er(M)S'"(PA^(M),^) = {Fh{(^),v) for all v e yn 

(25) "=' '"=' 

Y er(/x)6'"(9,wjv(/x)) = (G,(M),g) for ah q e Qh, 

m—l 

and the reduced basis equations (|20l) respectively into 

Ma A/6 

Y e:'{ti)A^iuNiti),v) + Y er(M)S'"(pAr(M),^) = {Fh{ti),v) for all v E yN 

(26) "=' "=' 



IN- 



Y ei''i(i)B'''iq,UN{(i)) = (G,(M),g) for all q G Qr 

m—l 

Once bases {^n}n=i ^^^ {in}n=i foi" the reduced basis spaces y^ and Qn respectively have been con- 
structed, the reduced system matrices 

[^"]m:=-4"(^L^£), l<k,£<2N; 
^'^^' [i3"]fe,, :=Z?"(|LCD- l<fc<2iV, l<i<N, 



are independent of the parameters /x, and can be precomputed and stored. This is the parameter- 
independent offline stage. In the online stage, for any fj. eT) we assemble and solve the 3-/V x 3iV system 
for the reduced velocity mjv G M^^^^ and reduced pressure p^ G R^^^ s.t. 



Ma \ /A/6 \ 

Y 0"(m)^" N^ + E 0r(M)s" p^ = j^(m) 

m—l / Vtti— 1 / 



(28) 

^ ' Mi 

Y^Tif^Wfu^^Giti). 

m—l 

In the case of pure Dirichlet boundary conditions we have 

Ma 

(29) F(m) = f''-Y 0a"(M)^'"«o 

m—l 

where the inhomogeneous Dirichlet condition has been treated by lifting the solution of the homogeneous 
equation into the space with proper boundary conditions. Therefore also the computation of the right- 
hand side splits affinely in terms of parametric dependence. The assembly and solution of the reduced 
system can be done independent (and thus very efhciently) of the dimensions of the finite element spaces 
for the velocity yh and pressure Qh respectively. While the offline stage is more expensive compared to 
solving the finite element problem, if the parametric PDE is evaluated for sufficiently many fj, (typically 
the cutoff point is between 100-500 PDE evaluations), the inexpensive online stage negates the preliminary 
costs involved in setting up the reduced basis matrices. 
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Figure 3. Result of the fluid-structure algorithm and a comparison of computational costs 



In our case the viscous transformation tensor vt obtained by free-form deformations does not satisfy the 
affine parameterization assumption (|24l) . This difficulty is handled by using the Empirical Interpolation 
Method (EIM) [l][9l[T9], which approximates the nonaffinely parameterized tensor with a suitable affinely 
parameterized one by replacing each term of the viscous transformation tensor with an approximate 
expansion 



Mi 



(30) 



Kk, = E^j^f(M)C(^)+£^'^(^,M), 



m— 1 



where the error terms are chosen to be under some tolerance, | |e*'-' (•, /i.) | |loo < e^j for all fj, & V. For this 
problem we use efJi'^' = le-5, which is satisfied with Mij < 22 for i,j = 1,2. The affine decomposition 
(f24| is then recovered. The tolerance should be chosen small enough so that the introduced additional 
error term does not dominate the error reduced basis approximation [121 119] . For an a posteriori error 
estimate of reduced basis approximation for the Stokes problem we refer at the moment to J18j . 



6. NUMERICAL EXPERIMENTS 

Because an analytical solution for the fluid-structure problem described in Sect. [2] is not obtainable, 
we verified our (finite element) computations by comparing them to results presented in flO] for the same 
test case and parameter values, except with a second-order wall law. For a flexible channel with length 
3 cm and half- width 0.5 cm we imposed a Poiseuille fiow profile at both the inflow and outflow 



(31) 



uo{x) = Vo 



1 - 4x2 




The numerical values for the physical parameters were chosen according to W so that the inflow velocity 
was vq = 30 cm/s and the blood viscosity i^ = 0.035 g/cm-s. The volume force on the fluid / was 
taken to be zero. The spring constant was K = 62.5 g/s^. From the initial guess fi — the least-squares 
optimization algorithm for the strong parametric coupling took 287 iteration steps (Stokes evaluation 



solution of least-squares problem for the coupling) until the stopping criteria ||/x 



fc+i 



^l- 



< le-5 



was fulfilled. The final value of the cost functional was J — 5e-8. The obtained displacement rj of the 



fluid-structure interface versus the assumed displacement fj are shown in Fig. 3(a) 



Once the full FEM solution was verified to be correct, we performed the model reduction of the Stokes 
problem. We consider two levels of model reduction. In the flrst level, the problem is parameterized 
with free- form deformations following Sect. [3], the transformation tensors are affinely decomposed with 
the empirical interpolation as in (pO)) , and the parametric decomposition (|25|) of the system matrices is 
precomputed. Then we can either: 

(1) Assemble and solve the full FEM Stokes-problem from ([25| using the precomputed system ma- 
trices for each fi. The only reduction is in terms of the parameterized geometry. We call this 
reduced FEM. 
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(2) Perform the full reduced basis offline stage and construct the reduced system matrices in (pS)) . 
then solve the reduced basis system for each /x. Reduction is both in terms of the parameterized 
geometry and the PDEs. This is the reduced basis approach. 

The computational costs are compared in Fig. |3(b)[ The costs are scaled so that the cost of one full 
FEM solution equals 1. The full FEM incurs no starting overhead, but scales poorly as the number of 
PDE solutions increases. The reduced FEM has moderate starting cost (empirical interpolation + matrix 
assembly) and scales fairly well. The full reduced basis approach is costly in terms of the offline stage, 
but solving the PDE is almost free, resulting in extremely good scalability when thousands of parametric 
PDE evaluations are desired. For the steady FSI problem we don't need thousands of PDE evaluations, 
so that reduced FEM would seem the best option in terms of total computational cost for this problem. 
For the unsteady case we anticipate a need for many more PDE solutions during time-stepping so that 
the added scalability of the full reduced basis becomes significant. 

7. CONCLUSIONS 

A model reduction technique for steady fluid-structure interaction problem of Stokes flow in a flexible 
2-d channel was proposed. A geometric reduction was performed using free-form deformations to reduce 
the free-boundary problem to a low-dimensional parameter space. Reduced basis methods were then used 
to reduced the complexity of the resulting parametric PDEs. 

A least-squares parametric coupling formulation between the fluid and structure was given. The 
approximate balance of stresses across the interface was formulated using the parameterized displacement 
of the interface and solved using nonlinear programming techniques. A low-dimensional parameterization 
with six parameters coming from a free-form deformation technique was enough to obtain approximate 
coupling. 

Computational costs between the full FEM, a geometrically reduced FEM, and the full reduced basis 
approach were compared. It was observed that the reduced FEM without the overhead of the reduced 
basis reduction was the computationally most attractive choice, but that for unsteady or nonlinear prob- 
lems there could be a need for the better scalability of reduced basis methods as the number of PDE 
evaluations increases. 
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